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ABSTRACT 

Equipartition arguments provide an easy way to find a characteristic scale for the mag- 
netic field from radio emission, by assuming the energy densities in cosmic rays and 
magnetic fields are the same. Yet most of the cosmic ray content in star-forming galax- 
ies is in protons, which are invisible in radio emission. Therefore, the argument needs 
assumptions about the proton spectrum, typically that of a constant proton/electron 
ratio. In some environments, particularly starburst galaxies, the reasoning behind 
these assumptions does not necessarily hold: secondary pionic positrons and electrons 
may be responsible for most of the radio emission, and strong energy losses can alter 
the proton/electron ratio. We derive an equipartition expression that should work in a 
hadronic loss-dominated environment like starburst galaxies. Surprisingly, despite the 
radically different assumptions from the classical equipartition formula, numerically 
the results for starburst magnetic fields arc similar. We explain this fortuitous coin- 
cidence using the energetics of secondary production and energy loss times. We show 
that these processes cause the proton/electron ratio to be ~ 100 for GHz-emitting 
electrons in starbursts. 

Key words: galaxies: magnetic fields - galaxies: starburst - radio continuum: general 



1 INTRODUCTION 

Magnetic fields are important in astrophysics, but in practice their strengths are very hard to directly determine. Their 
presence is inferred through synchrotron emission of cosmic ray (CR) electrons and positrons (e ) gyrating in magnetic fields 



Synchrotro n emission is detected from star-forming galaxies, demonstrating that they have magnetic fields (e.g., Condonlll992l : 
Beckll2005l ). However, the synchrotron emission only informs us of a combination of cosmic ray content and the magnetic field 
content. In a few cases for star- forming galaxies, a detected gamma-ray spectrum provides additional evidence of the CR energy 
density 1 Aharonian et al. 20061 : Abdo et al. 201ol and references therein), which when combined with p l ausible assumptions 
and modelling, allows us to co nstrain the magnetic fields in galaxies (e.g., de Cea del Pozo et al. 2009a : Crocker et al. 201ol . 
201fj : lLacki fc Thompson 2013]). Th e great majority of gala xies are undetected in gamma rays, though (some upper limits are 



given m 



Lenain fc Walterll201ll and Ackermann et al. 20121 ). Faraday measurements, which require a polarization signal, are 



a useful way of measuring ordered magnetic field strengths. Unfortunately, starbursts likely have highly turbule nt magnetic 
fields, so the expected polarizat ion is small. Anisotropy introduced in the t urbulence by shearing and comp ression (|Lainelll980l : 



Sokoloff et al1ll998l : lBeckll2012t may r esult in some polarization signal (c.f. iGreaves et al ] |2000l ; ljonesll2000h . but there still will 
be no Faraday signal (e.g., [Beck 2005) . Furthermore i n starbursts, Farada y depolarization can remove any polarization signal 
at typical observing frequencies (jReuter et al 19941 : Sokoloff et al. 19981 ). Faraday rotation measures also can be biased if 
there are relation ships between the magnetic field strength and density, or if anisotropic turbulence is present on the sightline 
( Beck et al.ll2003h . 



The equipa rtition and minimum energy arguments are common ways of finding a characteristic magnetic field strength B 
re.g., lBeckll200lh . For a given radio luminosity, if B is very small, then the CR energy density Ucr must be very large; likewise, 
if C^cr is very small, then B must be very high. In between, there is a single magnetic field strength where Ub = B 2 /(8iv) 
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is equal to Ucr'- this is the equipartition magnetic field strength (S cq ). Alternatively, for a given radio flux, the combined 
nonthermal energy density Ub + Ucr. has a min imum at a magnetic field strength -B m i n which is of the same order as (though 
generally distinct frorrQ) ^eq ( Burbidge 1956a V 

Throughout much of the Milky Way, equipartition between the magnetic fields and CRs holds (e.g., Niklas fcJ3eck] 



1997) 



Likewise, in other normal star-forming galaxies, equipartition likely holds to within a factor of ~ 10 in B (e.g.. lDuridll990l ). 
It is less clear if equipartition holds in starbursts, although any deviation would be interesting in its own right for th e 
propagation of CRs and the sources of the magnetic fields (e.g., Thompson et al. 20061 ; Lacki. Thompson. fc Quataert 2010l ). 
Equipartition has bee n proposed as a cau se of the observed co rrelation between the f ar-infrared and radio luminosities of 
star-forming galaxies IjNiklas fc Beckl ll997), On the other hand. iThompson et al. (2006) argued that equipartition (between 
magnetic fields and CRs) predicts magnetic field strengths too low to allow starbursts to lie on the correlation, since other 



radia tive losses are extremely fast and B must be large for thereto be any significant synchrotron emission (see also lLacki et al 



2010). However, the ease of equipartition methods - compared to methods involving gamma rays, detailed spectral modelling, 
or polarizati on methods - has led to their app lication in a variety of en vironments including normal galaxies, the Galactic 



Center (e.g.. lLaRosa et aL Il2005l; iFerrierj |2009| ). starburst galaxie s (e.g 



1 ^ J ' — - — — — — — '/ ' o ^7 

Persic fc Rephae li 2010; BeckJl2012j), and galaxies at high redshift 1 Murphvl 



Volk. Klein 
20091 ; IChakraborti et al 



Wielebinskilll989l ; beck et all 12005 
2012h . 

A problem with equipartition-style estimates is that most of the CR energy density is actually invisible in synchrotron 
emission: CR protons (and nuclei) are the dominant CR population at high energies. Only in the few cases where gamma- 
ray observations are available can the pro ton energy density be constrained directly ( Acero et al. 20091 ; Acciari et al. 20091 ; 
Lacki et alj|201ll ; IPersic fc Rephaelil [20121 ). Radio equipartition estimates therefore make assumptions about how to convert 
the observed CR electron population in some observed frequency range into the total CR proton energy density. The simplest 
and most common assumption is to assume that the CR proton and electron spectra are power laws in energy with the same 
spectral index, and that there is a single ratio k = 30 — 100 tha t sets the ratio for all galaxies. This is likely to be true for the 
injection spectra of primary CRs (iBelllll978l ; lschrickeiserj |2002), at least those ass ociated with star formatio n The spectrum 
of Milky Way CRs indicates the primary injection proton/electron is near ~ 100 (|Ginzburg fc Ptuskin|[l976l ). 

The problem with this assumption is that the CR spectrum can be complex, especially for CR . Roughly speaking, the 
steady-state CR spectrum N(E) is equal to the product of the injection spectrum Q(E) and the characteristic loss (cooling 
or escape) time t\ oss (E). When CR protons and electrons are governed by the same losses, such as diffusive escape in normal 
galaxies, the proton/electron ratio is preserved. In starburst galaxies, however, the CRs are likely to experience a varie ty of 
energy losses with different energy dependencies, which can alter the actual proton/electron ratio ( Beck fc Krausei 20051 ) . At 



high energies, e are cooled quickly by synchrotron and Inverse Compton (IC) losses, which increase s the proton/elec tron 



ratio with energy. This is seen in the Milky Way at energies > 10 GeV (see for example Figure 3.29 of ISchlickeiserll2002f l . In 
addition, ther e are ionization losses for both e ± and protons, bremsst rahlung losses of e ± , and winds, which can alter the CR 



spectrum (|Thompson et al.ll2006l : lMurphvlE)09l : lLacki etal1l2010h . 

A further complication is the possible presence of pionic secondary e in starburst galaxies. Secondary are generated 
when protons crash into ambient gas atoms, creating unstable pions that decay into gamma rays, neutrinos, and e . In the 
dense gases of starbursts, the amount of secondaries may be comparabl e to or even domin ant over the primary electrons (e.g., 
Torres! l20oi ; lRengaraia3l2005l ; IThompson. Quataert. fc Waxmanll2007l ; lLacki et allboioh . 

In this paper, we present equipartition and minimum energy formulae that should work in starburst environments. These 
take into account the presence of secondary and the strong energy losses of starburst galaxies. 



2 DERIVATION OF EQUIPARTITION AND MINIMUM ENERGY MAGNETIC FIELDS IN 
STARBURSTS 

2.1 Basic assumptions 

The spectrum N(E) of CRs is governed by the diffusion-loss equation; 



dt " ev ~ y 1 dE'" K ~'" K ~" t{E) ' ~ v " v "' ^ 

as given in Torres! 1 2004 ). Here Q{E) is the injection spectrum of CRs, b(E) = —dE/dt is the cooling rate of individual CRs, 
t(E) is the loss time to escape and losses that remove most of the CR's energy in one interaction (catastrophic losses), and 
D is the spatial diffusion constant. We can simplify the equation to the leaky box equation if we assume that the modelled 
region is homogeneous, so that all spatial terms drop out, and is in equilibrium (dN(E)/dt — 0). The leaky box equation is 
then 



1 When neither CR protons nor C R are cooled, the m inimum energy and equipartition magnetic field strengths are equal if CRs are 
injected with a spectral index of 3 iTBeck & Krausc 2 0051) . 
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(2) 



where escape out of the region is now considered a catastrophic loss. 

In the test particle approach to CR acceleration, CRs are injected with a momentum power law spectrum ( Belli 19781 ). 
To simplify matters, we assume that CRs are injected with an energy power law spectrum: 



, E \ ~ p 
Q(E) = Q — 



(3) 



where E is the total (rest plus kinetic) energy of the particle. The power law spectrum continues to a maximum energy 

7ma X mC 2 . 

At high ener gies, protons escape or lose energy through pion production. Since pion production can be modelled as a 
catastrophic loss (|Torresl l [200J), the proton spectrum is simply 



N P (E) = Q p (E)t p (E). 



(4) 



ine proton nietime m starbursts t p is thought to be dominated by advcction or possibly piomc losses m denser starbursts 
like Arp 220. The fact that starbursts are observed at TeV energies jAcero et al.ll2009l : lAcciari et al.ll2009h indicates hard 
gamma-ray spectra (r « 2.2), supporting the idea that t p is set by one of these energy- independent processes. The gamma-ray 
observations of M 82 and NGC 253 suggest that more CR proton power is advected away in their nuclear star burst winds 



than is lost to pion production, although the efficiency of pionic losses is still much greater than in the Milky Way IjLacki et al 



20 111 ). As we will see, though, it is convenient to scale t p with the pionic lifetime t^, because the secondary e injection rate 
is directly tied to the pionic loss time. We therefore parameterize the proton lifetime as: 

Tp = Fealty. (5) 

CR e ± in contrast are thought to largely be trappe d in starburst gala xies ( Volk 19891 ). although winds may be quick 
enough to remove before they cool in some starbursts ( Heesen et al. 201lh . Their losses are therefore continuous. The total 
e energy loss rate b(E) can also be written as a cooling time t e (E) = E/b(E), giving us a steady-state spectrum 

Q e (E)E Q e (E)t e (E) 



N e {E) = 



(p-l)b(E) 



To get the total CR proton energy density, we integrate its kinetic energy over all allowed energies: U p = 
Substituting in the solutions for the steady-state and injection proton spectra, we find 



Up = t^F cal Q (m p c ) T 
where we define a factor 



I _ -y 2 -P l-V" 
' /max . x /m; 



T 



(P/1,2) 



p-2 p-1 

(hl7 max + 1 - l/7max) (?? = 2) 



(6) 

N p (E)KdE. 
(7) 

(8) 



that depends on the shape of the injection spectrum. Unless otherwise noted in the paper, we assume 7 max = oo and 
parameterize p as 2.2 + Ap, so that T = 4.16 = 25/6 when Ap = 0. We can compare this to the luminosity of protons at 
energy E, which can be approximated as E 2 Q P (E). We find 

p-2 



tuf^i x e;q p (e) 



(9) 



2.2 Relating the secondary e to CR protons 



Calculating the secondary e ± spectrum invo l ves integrating the differential cross sections for pionic e ± 



production over the 

CR proton spectrum (e.g., Kamae et al. 20061 : Kelner. Aharonian. fc Bugavovll2006h . However, at sufficiently high 



e energies, 

the secondary e spectrum can be estimated using a S- function approximation. A fraction _F ca i of the CR proton energy goes 
into pion production, with the rest either escaping or lost to ionization or Coulomb cooling. Charged pions makes up ~ 2/3 
of the pions produced, with the remaining 1/3 going into neutral pions that mainly decay into gamma rays. Assuming equal 
ener gy going into the pion decay products, a pproximately 1/4 x 2/3 w 1/6 of the power in pions ends up in secondary e 
fc.f.. lBurbidgelll956bl : lLoeb fc Waxmanll2006h . The other half of the pionic luminosity (3/4 of the charged pion power) goes 
into neutrinos. For the 5-function approximation, we therefore have: 



2 E i 2 

E scc Qscc(E sec *) — — EpQ p (E p } 
b 



(10) 
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Table 1. Values of g for synchrotron emission spectra 



Loss process 


„ — d In t 
" — d In E 


g(p = 2.0) 


9(p = 2.2) 


g(p = 2.4) 


Synchrotron / Inverse Compton 


-1.0 


1.0 


1.00018 


1.0112 


Brcmsstrahlung 


0.0 


1.2343 


1.14299 


1.07964 


Ionization 


1.0 


3.00492 


2.26698 


1.83161 



Since the typical inelasticity per inelastic collision is 20 per cent, and each pion decays into four particles, the typical energy 



of a secondary e is roughly E s , 



0.05E P . 



Do secondaries in fact dominate the e spectra of starbursts? With equation[T0l we can estimate the primary to secondary 



ratio 
Q 
Q 



prim(^e) 6 \E e J 



(11) 



where we define Kq = Q p /Q pl i m as the number ratio of injected primary protons and primary electrons at high energies 
( > m p c 2 ). While the injection rate of protons vastly overwhelms that of primary electrons at a constant energy, the secondary 
are injected at much lower energy than their parent protons, where there is more power in p rimary electrons for p > 2. 



i_pr 

Thus the secondary e ± are 'diluted' with respect to the primary electrons fe.g.. lLacki et al.ll2010l ). For p = 2.2, we find 

Qscc [Ee) 
Qprim(E e 



6.9F ci 



(§) 



In terms of the fraction of the total e population in secondaries, / SC c = f 

-i 



cc/(Qsec H" Qprim) — [1 ~T" {Qscc /Qprim) 



fs, 



0.15F, 



al V 75 J 



(12) 



(13) 



for p = 2.2. 

The gamma-ray spectra of M 82 and NGC 253 indicate that F ca i 
imply secondary to primary ratios of ~ 1.4 — 3.4, or ,f s ec « 0.6 



0.2 — 0.5 IjLacki et al.ll201ll ). so our equations 
0.8. This is in line with detailed modelling of these galaxie s 
I Domingo-Santamaria fc Torreslliooil ; Ide Cea del Pozo. Torres, fc Rodriguez Marreroll2009bl ; iRephaeli. Arieli. fc Persidboifj ) . 



(14) 



2.3 Relating the radio emission to the e spectrum 

We also consider a 5-function approximation in calculating the synchrotron emission. In this approximation, we assume each 
e of energy E e radiates all of its synchrotron emission at a frequency: 

_ 3£feB(sina) _ 3E^eB 
Anmlc? 16m|c 5 ' 

with a magnetic field strength B and an electron electric charge e ( Rvbicki fc Lightman 19791 ). We have used the fact that the 
mean sine pitch angle (sin a) for an isotropic distribution of CR e is (sin a) = 7r/4. While physically speaking the S- function 
approximation is not correct since the synch rotron emission of an e ± is a broad continuum, it gives acceptably accurate results 
for a power law distribution of e ± energies ( Felten fc Morrison 1966J ). 

Since an e would radiate its energy in a synchrotron lifetime t sy nch, the volumetric synchrotron luminosity of a ho- 
mogeneous region is vt v = de/dlnv — (d(7e/dmj/)/t aync h = (E e dN e /dlnv)/t Byn ch, where v is the observed frequency. From 
ecmation ll4l dlni/ = 2dln_E e , so we have ve v w l/2(E e dN/dlnE e )/t Bync h = l/2(E^N e (E e ))/t Bync h- The steady-state spectrum 
of the electrons is given by equation [6] substituting that in, we get 



-E e Q e (E e 



gt e (E e 



(15) 



(p - l)t ay nch(£ e 

The factor g corrects our approximation for ve v to the exact result for a power law spectrum of ultrarelativistic e at the 
frequency vc ( Rvbicki fc Lightman 19791 ) . If the electron lifetime is dominated by a continuous energy loss process and scales 
as te(E e ) oc El , then the value of g is: 

27^3 



9 = 



-3(5+a-p)/2 7r (2+a-p)/2p 



^ 5+p~a ^ r ^ 3p-3a- 1 ^ p |- 19 + 3p-3a ^ p + 



(p-a+iy ~\ 4 r\ 12 r i 4 j r v 4 ji • ^ 16 ^ 

We compile the values of g for different loss processes and injection indices p in Table [T] but for bremsstrahlung, synchrotron, 
or IC cooling, g is very nearly 1, and it is exactly 1 when p — 2.0 and a — —1.0 (as in the case when synchrotron cooling sets 
the e ± lifetime). 

Now we can solve for the proton spectrum given a synchrotron volumetric luminosity. Combining equations [9] 1101 and 
1151 we have 
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12ff" 1 ' Ep 



p-2 



l)T^-y^-rtsyach(Se)/secI'ew. 



(17) 



The pitch-angle averaged synchrotron cooling time is 

67r(m e c 2 ) 2 



tsynch - ^£^2 . (18) 

where ot is the Thomson cross section ( Rvbicki fc Lightman 19791 ). Substituting it into eqn. 1171 and using eqn. [14] to put E e 
in terms of frequency and magnetic field, we find that 



U p = 



72imi e c 2 ( 20m £ 



gore V m p 



p-2 



V 3e 



V p_3)/2 B -^(p - i)T-^/ sec ^ 



te{E e 



(19) 



2.4 Which loss process is most important? 



A key ingredient that we need to derive the equipartition formula is the e cooling time. Losses from synchrotron, Inverse 
Compton (IC), bremsstrahlung, ionization, and advection all potentially contribute. The e lifetime from all of these processes 
is: 



[^synch ^IG 'brems "F ^ion F ^wind] 



synch 1" L IC f t broms Hon f "-wind] • (20) 

Our discussion here is similar to that of Thompson et al. 1 20061 ). Murphy (2009), and Lacki et al. 1 2010h . We consider 
M 82 as an example (the other nearby starburst NGC 253 has a similar environment), which ha s a magnetic field strength 
of ~ 200 iiG (from detailed mode l s of the radio spectr um; Dommgo-Santamarfa fc Torres! 20051: Persic. Rephaeli. fc Arieli 
2008; Ide Cea del Pozo et alj|2009bl : iRephaeli etldfaiol ), a mean gas densi ty of ~ 300 cm " ~ 3 JWeifi et alJboOll . assuming a 



scale height of 50 pc), and a radiation energy density of ~ 1000 eV cm ( Sanders et all 12003 . again assuming a starburst 



radius of 250 pc). Of course, none of these values are constant throughout the region, and it is not known how inhomogeneity 
affects CRs in starbursts, but using the average values seems to work for models fitting current radio and gamma-ray data 



I de Cea del Pozo et al.ll2009al ). After using eqn. [HI to convert between e energy to observed frequency, the synchrotron 



lifetime of GHz-emitting e is 



synch 



498 kyr 



B 



200 fiG 



-3/2 



(l GHz) 



-1/2 



(21) 



from equation ll8l For a radiation energy density !7 ra d the IC lifetime is tic = 3\/3/ (16(TT?7 ra d)-\/ m e ceB/v I Rvbicki fc Lightman 
19791 ): 



tic = 495 kyr 



B 



1/2 



(—Y 1/2 ( 

VlGHz/ V 



200 nG) VlGHzV VlOOOeVcm" 3 
The bremsstrahlung lifetime is only very weakly energy dependent, and at high e energies it is simply 

= 104 kyr ( n " 



(22) 



(23) 



300 cm- 3 , 

where we assume that nn = 10nHc in the interstellar medium ( Strong fc Moskalenko! 1998). The ionization lifetime is ti on = 



7/[2.7cgT(6.8 5 + 0.51n7)riffl, w here 7 = E e /(m e c 2 ) = *J 16m e cv / (3eB) is the Lorentz factor, and again assuming that 
nu = lOriHe ( Schlickeiser 20021 ) . If we evaluate the In 7 term as In 1000, since 1000 is the approximate Lorentz factor of 



GHz-emitting e in starbursts, then we find 



178 kyr 



B 



200 ^G 



-1/2 



1 GHz 



1/2 



300 cm' 



(24) 



Finally the advection time is roughly the time it takes to cross one scale height h at the wind spe ed « W md, with wind speeds 
of order a few hundred km s _1 observed in M 82 (e.g., iGrevel |2004 IWestmoquette et alJliooi) ! 2 ! Taking h = 50 pc and 
Vwind = 300 km s^ 1 , we find 



twind = 163 kyr 



50 pc 



1 300 : 



km s _1 



(25) 



2 We assume that both the gas and the synchrotron-emitting CRs have the same scale height h, since secondary e are created by 
interactions of CR protons with gas and the CR lifetimes are short in starbursts. 
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E„ [GeV] E„ [GeV] 

0.1 1 10 0.01 0.1 1 




0.001 0.01 0.1 1 10 100 1000 0.001 0.01 0.1 1 10 100 1000 

v [GHz] v [GHz] 

Figure 1. Comparison of the loss times from different processes in an M 82-like starburst (left) and in an Arp 220-like starburst 
(right). The total loss time t coo \ is the solid line, while bremsstrahlung is long-dashed, ionization is dotted, wind losses are long/short- 
dashed, and synchrotron and IC together are dash-dotted. On the top are the absolute loss times, while on the bottom is the ratio of 
each individual loss time to the total loss time. 



For our fiducial values of the physical parameters in M 82, we see that bremsstrahlung is the most important cooling 
process for 1 GHz emitting e , though ionization and winds are also important. However, the final e lifetime is t coo i « 39 kyr, 
shorter than even the bremsstrahlung losses. We can quantify the importance of other losses relative to bremsstrahlung with 
the ratio 

^brcms 



X = 



(26) 



For our fiducial M 82 values, X = 2.6 at 1 GHz; other losses are even less important. 

A comparison of how the individual loss times compare to each other at different frequencies is shown in Figure [1] (left) . 
In this figure, we include t he In -Tdependence in the ionization loss formula, and we use the full bremsstrahlung losses as given 



Strong fc Moskalenkol |l998)B Bremsstrahlung does dominate the losses of e that emit synchrotron between 300 MHz 



3 ( Downes fc Solomon 1998), and a radiation energy density equivalent to that of a 50 K blackbody (~ 

3 ). The loss times for these conditions are shown in the right panel of Figure [T] At 1 GHz, ionization now 

± 



and 6 GHz, with X ~ 3 in this range. At lower frequencies, ionization is the most important loss, while at higher frequencies 

synchrotron and IC are the most important losses. 

We also con sider the radio nuclei of Arp 220, using a magnetic field ~ 2 mG (e.g., Lacki et al, 2010l ). a density 
of ~ 10 4 cm 
30000 eV cm 

is the most important loss: this is because we see lower energy e ± with the higher magnetic field strength. Near 10 GHz 
though, bremsstrahlung again is the strongest loss, with X ~ 3. 

For other starburs ts, we generally exp ect advection, bremsstrahlung, or ionization losses to dominate: advection in low 
density starbursts (c.f.. lCrocker et al.ll201ll ). ionization in starbursts with high magnetic fields (likely in the densest starbursts 
like Arp 220) where we see low energy e , and bremsstrahlung in inter mediate density s tarbursts. At frequencies where 
bremsstrahlung dominates, X ~ 3 (compare with Figure 4, left panel of Lacki et al. 2010l ). In starbursts where advection 
is very strong, however, protons are likely to escape rather than produce pionic secondary e . Since advection is equally 
strong for protons and e ± , the proton to e lectro n ratio is just Ko when primaries dominate the e spectrum, and the revised 
equipartition formulae of lBeck fc Krausd (|2005l ) apply. 



2.5 The equipartition and minimum energy formulae 

For the purposes of this paper, we assume that the CR energy density entirely consists of protons (Ucr = U p ). Because there 
are several different energy regimes governed by different loss processes, the e spectrum is complex and likely has intrinsic 

3 The 7 < 100 for mula in equation C12 of IStrong fc Moskalenkol l|l998tl is erroneously multiplied by an extra m e c 2 , as can be seen in 
ISchlickeiserl {20021) , equation 4.4.17. We have corrected this when making Figure [l] 
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spectral curvature; accounting for all these loss processes would require knowledge of the density, radiation energy density, 
advection speed, diffusion constant, and magnetic field strength in the volume. More importantly, because there are always 
at least bremsstrahlung losses, and since tbrcms < iir, the steady-state proton/electron ratio will always be greater than the 
injection proton/electron ratio (c.f. equations [6] and [4]) . The injection electron/proton ratio is 



Q e (prim + sec) 1 
Q~v < 6 



V -Escc / 



-v 



(27) 



or < f 8 per cent for p = 2.0, < 11 per cent for p = 2.2, and < 7 per cent for p = 2.4. We therefore conclude that considering 
only protons leads to an error of less than ~ 20 per cent in U p . The other uncertainties are generally larger than this. 

Another possible worry of assuming that Ucn ~ U p is that the CRs include heavier nuclei, especially helium. In the Milky 
Way, though, the hydrogen/helium ratio at equal rigidities is ~ 8, and a t equal energies per nucleon, the hydrogen/helium 
ratio is ~ 24 ( Webber fc Lezniak 1974 ; Webber. Golden, fc Stephens! 1987 ). This ratio must be multiplied by the atomic mass 
A of helium when calculating the energy density, but still helium is only (1/24 — 1/8) x 4 ~ 1/6 — 1/2 of the CR energy 
density. Analysis of the CR flux at Eart h indicate tha t helium makes up about ~ 1/6 of the Galactic CR energy density, 
with another ~ 1/12 from heavier nuclei ( Webber 1998h . Furthermore, helium produces pionic secondary e ± just as hydrogen 
does: as far as collisions go, a helium nucleus at high energies is much like a collection of protons with the same energy per 
nucleon. Nucleons can shield each other in heavy nuclei, redu cing the pp cross section per nucleon as J 4 _1//4 for an atomic 
mass A jOrth fc Buffingtorj[l97l Istrong fc Moskalenkol l 19981 ). but for helium this just reduces the cross section by 30 per 
cent. Therefore, the hadronic starburst radio emission traces high energy CR helium nuclei as well as CR protons, and our 
derived U p basically includes the contribution from helium as well. 

The volumetric luminosity ve v can be found from either the total luminosity vL v or the intensity on a line of sight vl v . 
In the absence of absorption, the volumetric luminosity is simply: 



= vLu/V 



(28) 



where V is the volume of the radio wave emitting region. If the starburst is inho mogeneous, it is preferable to use the intensity 
instead of luminosity to calculate the equipartition field IjBeck fc Krausdl2005f ). The (unabsorbed) intensity on a line of sight 
is given by dl/d£ = e„/(47r), so that by assuming a constant e„, we can replace the volumetric luminosity with: 



where £ is the sightline through the radio emitting region. 



(29) 



2.5.1 Bremsstrahlung-scaled losses 

In general, the lifetime of the CR e is complex, depending on many loss processes. It is often convenient to scale the lifetime 
to the bremsstrahlung cooling lifetime. At GHz frequencies, bremsstrahlung is expected to be an important or even dominant 
cooling process for CR in starburst galaxies (section |2.4|I . Furthermore, both pionic and bremsstrahlung losses are mostly 
independent of energy, both scale with density, and both are very nearly of the same magnitude: 



U « 50 Myr(n/c 



from IMannheim fc Schlickeiserl l| 19941 ) for the pionic loss time, and (c.f. equation] 



31 Myr(n/cm 



for the bremsstrahlung lifetime. 

Equipartition formula - By setting B 2 /(87r) = CJJ P = C?7cr, the equipartition CR energy density is then 

4/(5+*) 



9m e c' 



goxc 

with a magnetic field of 

576tt 2 m e c 2 /20m 

-D, 



p - 2 f l6m e cu Y^y\ xu . ' 

( ^ J (P-1)Y- /soc^C 



2/(5+p) 



(30) 
(31) 

(32) 
(33) 



The ratio £ = Ub/Ucb. allows one to consider departures from exact equipartition. Here, X is the ratio of bremsstrahlung 
cooling time to £ C ool, as previously defined. Note that this assumes that X is constant in energy. However, X is likely to 
have some energy dependence, at least at energies far away from ~ 1 GeV. Since the frequency is a function of both energy 
and magnetic field, the magnetic field strength would enter implicitly into X; in order to solve self-consistently for B, the 
energy dependence of X must be included when solving eqn. 1191 In practice, because of the weak dependence of B eq on X 
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Table 2. Ratio of equipartition and minimum energy magnetic fields 



Formula 



B c . 



B„ 



B„ 



-(p = 2.0) 



B„ 



-(p = 2.2) 



B G , 



B„ 



-(p = 2.4) 



Reference 



Classical formula 
Revised formula 



[4/3] 2 / 7 
[4/0+ l)] 2/(p+5) 



1.09 
1.09 



1.06 



1.04 



Beck fc Krause (2005) 
Beck fc Krause (20051 



Hadronic formulae 



Bremsstrahlung-scaled losses 
Ionization-dominant losses 
Synchrotron-dominant losses 
IC-dominant losses 



[4/0 + l)] 2/(p+5) 
[ 4 /p]2/(p+4) 

[4/(p-2)] 2 /(P+ 2 ) 
[4/(p + 2)] 2 /(P+ 6 ) 



1.09 
1.26 
oo a 
1.0 



1.06 
1.22 
4.2 
0.99 



1.04 
1.17 
2.8 
0.98 



§ [2311 (Eqns. [33lfc[36t 
§ l2~ST2l (Eqns. I42lfc l45)l 
§ [2~Q1 (Eqns. l49lfc[52t 
§ [23^1 (Eqns. [58lfc[6ll 



a : As discussed in Section l2.5.3l in practice this limit can never be attained, because of breaks in the CR spectrum, non steady-state 
effects, and other radiative losses that become important when B is low. 



(B (x X 2/7/2 ),the effects of energy dependence in X on the estimated equipartition magnetic field will be small. Furthermore, 
below the ~ GeV energies probed by GHz observations, dX/dE < because ionization losses are quicker at low energies, and 
above ~ GeV, dX/dE > because Inverse Compton and synchrotron losses are quicker at high energies, so we roughly expect 
dX/dE £s (as shown in Figure [Q. 

Suppose that the starburst is a homogeneous disc, and the synchrotron-emitting particles fill a fraction f^ ch of that 
disk. If it has radius R and midplane-to-edge scale height h, then the emitting volume is V = 27ri? 2 /i/fljf\ The observed flux 
density S„ is L^/(4tvD 2 ), where D is the distance. Scaling to values typical of nearby starbursts like M 82 and NGC 253: 



U C r 



106 eV cm' 



(l GHz) 



2.4+2Ap 
2 + Ap 



23.4 



1 + 



.X 



19/ rsynch 
' /fill 



3.2 + Ap 

" 7.2+Ap x 



( R 



65.4 



" G (nk) 



1.2 + Ap 
7.2 + Ap 



23.4 



1.2 J 



fsccX 

rsynch 
/fill 



^250 pc 



( R 



It 



100 pc 



D 



( T 



3.5 Mpcy ^25/6 



^250 pc 



h 



100 pc 



D 



3.5 Mpc 



(25/5) 



g 
1.14 



9 
1.14 



7.2+Ap 



(34) 



7.2 + Ap 



(35) 



Minimum energy formula - The minimum energy magnetic field is defined as the magnetic field which minimizes U b + Ucr 
for the observed synchrotron luminosity. It is found by setting d[B 2 /(87r) + U p ]/dB = 0. Numerically, the minimum energy 
estimate is very nearly the same as the equipartition estimate: 



B n 



O + i) 

61.3 fiG 



144tt m e c 
g<JTC 



20m £ 



16m e ci/\( p - 3 )/ 2 



3e 



(P- 



f-brems 



2/Cp+S) 



1 GHz 



1.2 + Ap 
7.2 + Ap 



24.2 



Ap 



(-i)( 



X 



3. 2 J f syn* 



R 



250 pc 



100 pc 



D 



I- 



T 



3.5 Mpc J ^25/6 



(i 



(36) 



7.2+Ap 
(37) 



The ratio B cq /Bmi n is [4/(p + 1)] 2 /(p+ 5 ) i the same as in the revised equipartition formula of Beck fc Krause (2005) (see 
Table [2}. When bremsstrahlung losses dominate, the proton to electron ratio is constant with energy, so the proton/electron 
ratio has the same behavior as if there are no losses. 



2.5.2 Ionization-dominant losses 

At low energies, ionization (or Coulomb losses) will always dominate the energy losses of CR . For M 82, ionization losses 
should become most important for e radiating below ~ 300 MHz, while for Arp 220, ionization losses are most important 
all the way up to ~ 3 GHz (Figure [TJ. This assumption is therefore most appropriate for observations with low frequency 
telescopes like LOFAR when they observe starbursts. While more complicated to scale from pionic losses, because the ionization 
losses depend on electron energy, at sufficiently low energies they will be the only loss, making the calculation of the electron 
spectrum relatively simple. 
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tion — C 



Ignoring logarithmic terms in energy, the ionization loss time for CR e can be parameterized as: 
B \- 1 ' 2 /i/\V2 



V-Bo/ \VqJ \cm J / 



(38) 



where ii on ,o = 53 Myr for Bo = 200 /iG and vq = 1 GHz (compare equation I24[) . Likewise, the pionic loss time can be 
parametrized as: 



Vcm J / 



This lets us cancel out the density dependence in the ionization and pionic loss times. 
From equation 1191 we have 

p-2 



U p = 



727rm e c^ 
ga T c 



\ m p J V 3e / ^ t ion , V-Bo/ V^o/ 



SCC 



Equipartition formula ~ Solving for the equipartition magnetic field B eq = V^C^cr, we find 



U C R 

and 



-B eq = 



72nm e c 
gaTC 



(87T) 



-p/4 



\ 37 



Mon,0 Wo J 



-p/4 



4/(p+4) 



5767r 2 m e c 2 / 20m 



garc 



p-2 



1 2/(p+4) 



(39) 



(40) 



(41) 



(42) 



For typical starburst values, we have 



C/cr 



33.6 eV cm 



0.4+2Ap 
V \ 6.2 + Ap 2.2 + Ap 

£ 6.2 + Ap X 



(l GHz) 



Ap 



31.2 
and 

B cq = 36.8 jUG 



1.2/ /=y nch ^ Jy y ^250 pc 



h 



D 



100 pc J V3.5Mpc/ V 25/6/ V 2. 27 



/ T 



6.2 + Ap 



(43) 



0.2 + Ap 
f \ 6.2 + Ap 



1 GHz 



31.2 



Ap 



i + ( f^ x (!k\( R 
1.2 J /^ nch v J yy V 250 p c 



£> 



100 pcy l3.5Mpcy ^25/6yV2.27 



6.2 + Ap 



(44) 



Minimum energy formula - As with the bremsstrahlung loss formula, we solve for B m i n by setting d[B 2 / (fin) + Up\/ dB = 0. 
We find 

2/(p+4) 

(45) 



1447r 2 m e c 2 / 20m 



go~TC 



V 3e / ti on ,o Vfo/ 



SCC ^ C I> 



Numerically, this comes out to 
B min = 30.4 fj,G 



0.2 + Ap 
1/ \ 6.2 + Ap 



1 GHz 



34 



.Ap / Ap W Ap\ / SCC X / A 

V 1.2 y V 2.2 J I Jy 



R 



/mi V J y/ V 250 pc / V 100 pc / V 3 - 5 Mpc 



D 



(jl.) f^r 1 

^25/6 J 12.277 



6.2 + Ap 

(46) 



2.5.3 Synchrotron- dominant losses 

At high energies, synchrotron and Inverse Compton losses always dominate the cooling of CR e . For M 82, this should 
happen for e emitting above ~ 5 GHz, and for Arp 220, the transition is for e emitting above ~ 15 GHz (Figure [TJ. 
In order for the observed infrared-radio co rrelation to hold for starbursts, it is thought that Ub ~ U Ta d in starbursts l|Volkl 
19891 : ICondon et al.lll99ll : lLacki et alJboioh . Therefore for high freq uency radio emission , it ma y be more convenient to scale 



to the synchrotron lifetime. This case was previously considered in lPfrommer fc Enfilirj (|2004T ). but we rederive it with our 
more approximate approach to compare with the other formulae. Synchrotron losses do not depend on density, so the density 
dependence in t w is not canceled out. 
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As we can see from equation 1171 t coo i = t sync h cancels the t sync h in the CR e ± emissivity. After converting the electron 
energy E e w E p /20 (where a pionic e ± typically has 1/20 the energy of its primary proton; section T2.2p into synchrotron 
frequency using eqn. 1141 we are left with: 

p-2 



U p = Ug- 1 ^ 



x / 16rn e cv 



(p-2)/2 



13 



-(p-2)/2 



/ 20m £ 



(p - 1)T^/ S 



sec v c (y 



3e / y m p 

Equipartition estimate - Setting B 2 /(8ir) = £(7p = £Ucn, we find 



C/cr 



12<7 



V 3e / y m p J 



(p-2)/4 



4/(p+2) 



The equipartition magnetic field strength is 



96ng 



! /16m e c//\(f- 2 )/ 2 /20m e \ P 2 



l)Tf w / S ec^e^C 



2/(p+2) 



(47) 



(48) 



(49) 



Substituting in the previously used values, we find 



U C r 



5.45 eV cm 



(l GHz) 



4.4 + 2Ap 

4.2 + Ap 0.2 + Ap 

£ 4.2 + Ap X 



,Ap 



49.2 

for the energy density and 



Ap\ f BCC X 
l 2 ) ff ch 



250 pc 



h 



100 pc 



3.5 Mpc 



( nH V 

V300 cm" 3 / 



/ T 



V 25/6 



4.2 + Ap 

(50) 



14.8 



" G (nk) 



2.2 + Ap 
4.2 + Ap 



49.2 



j-synch 
fill 



Jy / ^250 pc 



ft 



100 pc 



3.5 Mpc 



300 cm ~ 3 



— 1 / 'Y* 



^25/6 



1.0 



4.2 + Ap 

(51) 



for the equipartition magnetic field strength. 

The resulting energy densities and magnetic fields are much smaller than the bremsstrahlung estimates. This is because 
the synchrotron cooling alone is relatively slow compared to pionic cooling for the parameters used. Thus, the pionic e 
accumulate longer, leading to a smaller proton/electron ratio. With less 'invisible' CR proton content compared to CR e 
content, the energy density is smaller. 

Minimum energy formula - The minimum energy magnetic field strength with these assumptions is 

_ 2 -1 2/(p+2) 

247rg (p-2)^ — — — J l^ - ) (p-l)T^/ scc ^ (52) 



With our usual parameters, we find that B m i n is much smaller than B c , 



3.56 



2.2 + Ap 
4.2 + Ap 



100 



D 



3.5 Mpc 



f ng V 

V300 cm" 3 / 



h 



+ fsccX (^L 

\ 1.2/ V 0.2/ /=y nch y^Jyy ^250 pcy V 100 P c 

2 

— )( JL Y 1 ' 

25/6 I Vl.O/ 



(53) 



The reason that B m \ n <C -B eq (Table [2]) is that the derived U p depends very weakly on B. When all of the electron power 
is going into synchrotron, we have ve v = E 2 N e (E e ) It^^Yi = -EeQe(-Ee)£synch/tsynch = E 2 Q(E e ). Even as B gets arbitrarily 
small, the assumption that e ± cool only by synchrotron implies that, in steady-state, the synchrotron luminosity is exactly 
equal to the injected power. The only way that B affects the synchrotron flux is the dependence of frequency on magnetic 
field strength: varying the magnetic field strength probes different parts of the e spectrum, where there are different amounts 
of power injected if p ^ 2. This is a weak effect. Therefore, since the derived U p is nearly independent of B, the nonthermal 
energy density can usually be decreased by simply decreasing Ub with little effect on U p . When p = 2.0 exactly, then i? m j n 
goes to 0. With ap = 2 spectrum, there is equal power injected at all energies, so the spectral effect disappears. In practice, 
this limit is not attained because (1) there is a high energy cutoff to the CR e spectrum, and the frequency it is observed 
at depends on B; (2) steady-state may not be attained if synchrotron is the only radiative loss; and (3) as B decreases, the 
other losses become more important; if nothing else, when B < 3.5(1 + z) 2 fiG, Inverse Compton losses from the CMB exceed 
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synchrotron losses, so this estimate is no longer valid. Finally, if synchrotron really were the only losses, then as £ sync h — > oo, 
U e becomes bigger than U p , invalidating our assumption that Uck ~ U p . In practice, this cannot happen for hadronic 
because bremsstrahl ung losses would always int ervene (as discussed in the beginning of the section). 

For comparison, IPfrommer fc Enfi lin (2004) found a minimum energy magnetic field strength for pionic e ± that can be 
written as 

1 2/(p+2) 



Bi 



2567r 2 (p-2)(p + 2) B(^±, *?)r(*J5) 



9VMp + io/3) r(^)r(^±iH)r(£±^) V m 



16m e 



p-2 



1 2nm e cv 
V 3e 



(p-2)/2 



(54) 



when synchrotron losses dominate and p / 2, with a vv as the cross section for proton-proton pionic collisions (given in their 
equation 2.18), B is the beta function, and F is the gamma function. This version of the estimate assumes there are no primary 
electrons (/ sec = 1). Plugging our fiducial parameters into this equation instead of equati on 1521 gives us 3.7 /iG instead of 
3.6 fj,G in ean. 1531 The result from our approach is therefore confirmed bv lPfrommer fc Enfilinl (2004). 



2.5.4 IC-dominant losses 



The other dominant loss at high energies is Inverse Compt on (IC). In starburs t galaxies, IC losses are expected to be quick 
because of intense FIR radiation from star-formin g regions (|Condon et al.lll991h ; the only question is whether the other losses 
are stronger (see the discussion in [2^llVolk|[l989T ). The IC loss time has a convenient scaling with the synchrotron loss time: 



tlC = isynch(t^s/f7rad) 

where [/ ra d is the radiation energy density. We therefore have 



(55) 



_i (lQm e cv 



(p-2)/2 



B 



-(p-2)/2 



( 20m, 



p-2 



(p — l)T£,r/sec-77^f tv 



(56) 



These equipartition and minimum energy estimates require both a density and a radiation field. 

Equipartition estimate - By setting the ratio of magnetic energy density to CR proton energy density as we find: 

p-2 1 4 /(f+ 6 ' 



£/cq = 

and 

Boa = 



Ug- 1 (8tt) 



_ 2 )/4 ( 16Tn e cv \ ( p ~ 2 >/ 2 ( 20m, 
V 3e 



7Rfi 2 -1 ^ 16m e c^ V"- 2 )/ 2 ( 20m e V \ , m . „ , 
V — 3e — ) \~m — J ~~ l)Tt w / 80 cl/tadfevC 



2/(p+6) 



with the accompanying scalings of: 



U C r 



69.3 eV cm" 



D 



(l GHz) 



4.4 + 2Ap 

8.2+Ap 4.2 + Ap 

t 8.2 + Ap 



3.5 Mpc 



n n 



300 cm 



' T ^ (A)"' 



1000 eV cm~ 3 7 ^25/6 



Jy J \ 250 pc 

4 

8.2+Ap 



100 pc 



and 



B eq = 52.8 



2.2 + Ap 
8.2 + Ap 



— V( 

3.5 Mpc y 1 300 cm" 



26 l Ap (l I ApN | C fsccX ( — 



/ J? 



Jy y V 250 P c 



100 pc 



/ VlOOO eV cm-V ^25/6 / V 1.0/ 



8.2 + Ap 



Minimum energy estimate - Using eqn.UHlto find the minimum nonthermal energy density, we calculate 



-^min 



2/(p+6) 



(57) 



(58) 



(59) 



(60) 



(61) 



Bmin = 53.5 /lG 



2.2 + Ap 
V \ 8.2+Ap 



1 GHz 



Jy 7 V 250 P c 



100 pc 
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D 



\ 2 

)G 



mi 



3.5 Mpc I V 300 cm 



Ur,d \ ( J_\ (J_ 

1000 eVcm- 3 7 I 25/6 J 1 1.0 



8.2 + Ap 



(62) 



for p = 2.2 + Ap. For p near 2, the minimum energy and equipartition magnetic field strengths are nearly identical in this 
case (Table El). 



Pfrommer fc Enfilinl [|2004h also derived the minimum energy estimate for IC-loss dominated pionic e in a more sophis- 



ticated manner, finding: 



R P04 



2048^ 3 (p + 2) 2 B(£^,^E)r(£±5) 

9^( P + io/3) r(^£±2)r(^±i2)r(2±^) 



16m e 



p-a 



27rm e ci/\( p - 2 )/ 2 z>e„£/ rad 



V 3e 



2/0+6) 



(63) 



Here, we have substituted the general radiation energy density [7 ra d for the CMB energy density. Plugging our fiducial starburst 
parameters i nto this equation gives us 5 5 pG as the coefficient for equation 1621 Like the synchrotron- loss dominant formula, 
the results of iPfrommer fc Enfilinl (120041 ) are in accord with our own. 



3 APPLYING THE EQUIPARTITION FORMULA 

3.1 How much of the radio flux is diffuse and nonthermal? 

So far, we have been assuming throughout the paper that all of a starburst's GHz radio flux is diffuse synchrotron emission 
from CR e^ 1 in the interstellar medium. We also assumed that the synchrotron flux is transmitted freely to Earth. Neither is 
completely true, although it turns out these are often reasonable approximations. 

At high frequencies, thermal free-free emission, which does not trace magnetic fields or CRs, becomes increasingly im- 
portant, since it has a flat v~ ' 1 spectrum. However, the total amount of free-free emission is limited by the number of 
ionizing photons generated by the starburst. At 1.4 G Hz, the therma l fraction is at most ~ 1/8 for starbursts that lie on 
the observed infrared-radio correlation, as estimated by I Condon (1992): if many ionizing photons escape the starburst or are 
destroyed by dust do not contribute, the thermal fraction might be sma ller. A few very you ng starbursts appear to have no 
synchrotron flux, perhaps because supernovae have not gone off in them ( Roussel et al. 20031 ). but most star-form ing galaxies 
and starbursts have radio spectral indi c es <C —0.1, indicating nonthermal em i ssion dominates the GHz flux (e.g.. iKlein et al 



19881 : ICondorj[l992l iNiklas et all 1 19971 : Iciemens et al.ll2008l : llbar et al.l bold ; I Williams fc Bowerlboid). Radio spectrum fits 
to normal spiral galaxies indicate their typical 1.4 GHz thermal fraction is ~ 10% (jNiklas. Klein, fc Wielebin ski 1997). For 
starburst galaxies, the spectral fits are more difficult, because the t ransition between the d ifferent cooling processes at different 
frequencies cause the intrinsic synchrotron spectrum to steepe n dThompson et aL 20061 ). whereas free- free emission flattens 
the spectrum; disentangling the effects of the two is difficult (Condon 1992 ). Model fits to starburst radio spectra so fa r 
generally indicate GHz thermal fractions of a few percent (e.g., Klein. Wielebinski. fc Morsil 19881 : Williams fc Bower 2010l ). 



We adopt /therm = [9(i^/GHz) _0,6 + which is appropriate if the nonthermal synchrotron spectrum has a v~ ' 7 spectrum 
and the thermal fraction at 1 GHz is 10%. 

The same ionized matter that generates free-free emission must also have free-free absorption on some level, but the 
effects of absorption are harder to calculate because it depends on the geometry of the ionized gas. The simplest approach to 
estimate the free-free absorption from the radio spectrum is with uniform slab models, in which the synchrotron-emitting CRs 
and ionized 10 4 K gas are assumed to be homogeneous and fully mixed. With uniform slab models, the free-free turnovers are 
generally expected to be in the range of a few hundred MHz for M 82 and NGC 253's starburst cores, with optical depths of 
order ~ 0.1 at 1 GHz (|Carillilll996l : lwilliams fc Bowerllioiol : Udebahr et al.ll2012l). Free-free absorption may be important for 



frequ encies up to a few GHz in Ultraluminous Infrared Galaxies (ULIRGs, Condon et al. 1991 : Torres! 2004 : Clemens et al, 
2008). In either case, free- free absorption might significantly suppress the observed synchrotron luminosity at frequencies 
where ionization cooling dominates 

Yet the ionized gas is not likely to be uniformly distributed in the starburst. Much of the ionized mass instead is likely 
concentrated discrete HII regions, which are dense but have small filling factor. Since not all sightlines through the starburst 
necessarily pass thro ugh an H II r egion, it is possible some synchrotron flux is transmitted even as the frequency descends deep 
into the MHz range jLackilboij ). Radio recombination line observations are often interpreted with models of H II regions. 
These studies find H II regions of a wide range of densities in starbursts, but generally fall into the categories of relatively 
low density (~ 1000 cm~ 3 ) regions wi th relatively high covering fractions and relatively high density (~ 10 5 cm" 3 ) region s 



with relatively low covering fractions ( Anantharamaiah et al. 19931 : Rodriguez-Rico et al. 2004 : Rodriguez- Rico et al. 20061 ) . 
In Arp 220, the high density HII regions are opaque even past 10 GHz, but they cover onl y a fraction of a percent of the 
starburst; a more diffuse component of ionized gas has an optical depth of ~ 10% at 8 GHz ( Anantharamaiah et al. 200d ) . 

Given the difficulty of accurate calculations of free-free absorption, the relatively weak dependence of the equipartition 
and minimum energy estimates on the flux, the uncertainty in other factors like emitting volume (both because of projection 
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Table 3. Updated equipartition estimates for starbursts 
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Bremsstrahlung-scaled" 
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(Mpc) 


(pc) 
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M 82 


3.6 


250 


50 


1.0 


8.94 


940 


190 


180 


160 


240 


220 


(1) 


I\ (_/(_. 253 Oore 


3.5 


150 


50 


1.0 


3.0 


880 


190 


180 


160 


230 


220 


I 1 ) 


IN GO 4945 


3.7 


540 


50 


1.4 


4.2 


300 


110 


100 


89 


130 


130 


( 2 ) 


NGC 1068 Starburst 


13.7 


3000 


50 


1.4 


1 


86 


59 


55 


47 


72 


68 


(3) 


IC 342 


4.4 


710 


50 


1.4 


2.25 


190 


87 


82 


71 


110 


100 


(1) 


NGC 2146 


12.6 


520 


50 


1.4 


1.09 


570 


150 


110 


130 
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180 


(1) 


NGC 3690 


42.2 


2460 


50 


1.4 


0.66 


300 


110 


100 


89 


130 


130 


CI) 


NGC 1808 


14.2 


1030 


50 


1.4 


0.52 


200 


91 


85 


73 


110 


100 


(4) 


NGC 3079 


15.9 


190 


50 


1.4 


0.85 


2000 


280 


270 


240 


350 


330 


(4) 


Arp 220 West 


75 


100 


50 


8.4 


0.072 


9300 


610 


580 


500 


750 


710 


(5) 


Arp 220 East 


75 


100 


50 


8.4 


0.061 


8500 


590 


550 


480 


720 


670 


(5) 


Arp 193 


99 


160 


50 


8.1 


0.035 


5000 


450 


420 


360 


550 


520 


(4) 


Mrk 273 


160 


160 


50 


8.1 


0.044 


9800 


630 


590 


510 


770 


720 


(4) 



Columns: D is the distance to the starburst; R is the assumed radius of the starburst disc and h is the assumed midplane-to-edge 
scale height of the starburst disc, for a total volume of V = 2nR 2 h; S„ is the total observed flux density at frequency v. We take 
the thermal fraction of the radio emission to be /therm = [9(^/GHz)~ !i + 1] — 1 . 

a : Equipartition energy densities, equipartition magnetic field strengths, and minimum energy magnetic field strengths calculated 
assuming radio emission is hadronic and with lifetimes scaled to the bremsstrahlung loss time, from equations 1321 1351 and 1361 
We assume X = 3, / scc = 1, f^f = 1, £ = 1, and p = 2.2. 

: Classical minimum energy estimate of magnetic field strength from lLongairl j2010t ). which assumes Ko = 100, Vmin = v, and 
a = 0.75. 

c : Revised equipartition and minimum energy magnetic field strengths from iBeck fc Krausd l|2005h . which assumes p = 2.2 
(synchrotron spectral index a = 0.6) and Ko = 100. 

d : References - (1): M 82 and NGC 253 r adio flux densities from IWilliams fc Bowerl l|2010l l: (2): NGC 4945 radio flux 
density and size from IStrickl and et al.l fl2004h ; (3): Radio flux density and approximate radius of NGC 1068 st a rburst from 
IWvnn- Williams. Becklin, &c Scovill el [119851 ); (4): Radio flux densities and radii compiled fro m iTh ompso n et al. I j2006h ; (5): 
Distance and sizes of Arp 220's radio nuclei from lSakamoto et alj {§008); radio flux density from IThompson et al. I (l2006fT 



effects and because of the unknown value of the CR filling factor), and the fact that free-free absorption optical depth is 
expected to be significantly less than 1 at the frequencies we consider (1 - 1.4 GHz for most starbursts, 8.4 GHz for ULIRGs), 
we ignore it when estimating the magnetic fields and CR energy densities below. However, its effects should be considered 
carefully when working at low frequencies with the ionization-dominant losses formulae. 

Finally, not all of the synchrotron flux comes from the CRs in the starburst interstellar medium. First, some of the radio 
emission actually comes from individual radio sources like supernovae remnants. The fraction of the flux from these individua l 
sources is only a few percent of the total starburst flux, though, so we ignore it (|Lisenfeld fc Volklbood : lLonsdale et al.l l2006). 
Second, not all of the flux comes from the starbursting region itself, where the CRs are accelerated or created through pion 
productions, and where the magnetic and CR energy densities are though to be highest. Some comes from the surrounding 
host galaxy, and additional emission can arise from a "halo" region thought to arise when CR are advected away from 
the galaxy. In the case of NGC 2 53, about half of the radio flux in fact com es from these regions rather than the starburst 
core, so we use the core flux only jHeesen et al.ll2009l ; Iwilliams fc Bowerlboid ). Likewise, in NGC 1068, most of the radi o flux 
comes from an active nucleus and its jets, so we use the estimated flux of the starburst alone (|Wvnn- Williams et al.l ll985). We 
assume in the other cases that all of the synchrotron flux comes from the starburst core, but this can be tested with resolved 
observations. This is known to be a re asonable approximati on for M 82, where the starburst core does actually dominate the 
total flux at frequencies above 1 GHz l Adebahr et al. 20121 ). 



3.2 New equipartition estimates of B in starbursts 

We present new equipartition estimates of the magnetic field strength in selected starburst galaxies in Table|3] These estimates 
use the bremsstrahlung-dominated formula with X — 3, and also assume that p — 2.2 and = / scc = £ = 1. For comparison, 
we also give the results for the classical and revised equipartition formula. Note that we use 8.4 GHz radio data for the ULIRGs 
Arp 220, Arp 193, and Mrk 273; at these higher frequencies, bremsstrahlung still should be the most important loss (Figure[T] 
right panel). 

The equipartition magnetic field strengths range from 60-600 jj,G. For M 82 and NGC 253, t he estimated magnetic 



fields a re 200 nG, which i s enti rely in line with broad band modelling of their nonthermal spectra (e.g.. lde Cea del Pozo et al 



2009a ). Klein et al. 1 19881 ) and Adebahr et al. 1 20121 ) found smaller magnetic field strengths of 50-100 /xG, but they assumed 
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the emission was spread over a larger volume so that the energy densities were smaller. The magnetic field strength estimates 
for the distant ULIRGs such as Arp 220 are of order half a milliGauss, which is several times smal ler than the few milliGauss 
expected from the linearit y of the FIR- radio correlation ( Condon et al, 1991 ; Lacki et al. 2010T I and detailed modeling of 
Arp 220's CR population (jTorres! 120041 ). The derived magnetic field strengths depend weakly on £ = Ub/Ucr, though, 
with B oc £ 2 /( 5+p ' and [/cr oc £ _ ( 1 +w/( 5 +i>) for the bremsstrahlung-loss formulae. At a given radio flux, and if no other 
assumptions are changed, magnetic field strengths that are 3 times stronger than equipartition require £ w 50 and cosmic ray 
energy densities 6 times smaller than the values derived for £ = 1. 

However, the most interesting feature of the new equipartition estimates evident in Table [3] is that they are still of the 
same order as the classical and revised equipartition and minimum-energy estimates. 



K 



I E p \ p ~ 2 
V E scc J 



3.3 Why do previous equipartition formulae seem to work? 

If previous equipartition formulae do not apply to starbursts, because of starbursts' strong electron cooling and the presence 
of secondary e , why do these formulae give similar results to the corrected formula? 

The answer is that the steady-state proton/electron ratio at ~ GeV energies is not far from the canonical ~ 100 in 
starbursts, although the reasons for this factor are more subtle than usually realized. The steady-state proton/electron 
number ratio at any given energy is 

NpjE) = 6UX(p - l) /sec ( E p y~ 2 

rerr 

In a starburst like M 82 with F cai = 0.4, p = 2.2, and X = 3 (and therefore / scc = 0.73), we find K = 120. Even for F ca i = 1, 
K = 57. Therefore, an estimate of the proton energy density that assumes a constant proton/electron ratio of K = 100 is 
actually correct to order of magnitude for 1 GHz observations. 
The reasons the assumption works are that: 

• While synchrotron and IC strongly suppress the proton/electron ratio at high energy, and while ionization losses strongly 
suppress the proton/electron ratio at low energy, bremsstrahlung-dominated losses roughly preserve the proton/electron ratio, 
since ^Tr/farems ~ 1-6. Although bremsstrahlung is more important at high density, reducing the steady-state electron density 
( Beck fc Krausell2005h . -he pionic losses also scale with density, reducing the steady-state proton density. 



• At 1 GHz, we are observing at intermediate energies in starbursts, a sweet spot where none of the other losses greatly 
overpowers bremsstrahlung. But while bremsstrahlung may be the strongest individual loss, the other losses combined are 
more important than bremsstrahlung (X « 3). As seen in Figure [T] the e radiative lifetime reaches a maximum for energies 
observed at GHz frequencies. So the combined effect of the losses is to suppress the proton/electron ratio, but only by a factor 
of ~ 5. 

• On the other hand, the proton/electron injection ratio is raised by a factor ~ 5 due to the presence of pionic secondaries. 
The 'dilution' effect arising because secondary e are injected at much lower energy than the primary protons helps ensure 
that the secondaries do not swamp the electron population. 

This is related to the 'high-£ s conspiracy' posited by Lacki et al. 20ld to explain the FIR-radio correlation from normal 



galaxies to starbursts: radio emission is enhanced by secondaries, but is suppressed by non-synchrotron losses including 
bremsstrahlung. 

A consequence is that if we move to higher or lower electron energy (or observing frequency), X will increase beyond 
~ 3. The coincidence no longer works, and the proton/electron ratio is increased by either ionization losses at low frequency 
or synchrotron and IC losses at high frequency. At sufficiently low energies ( < 100 MeV; observed < 40 MHz in the M 82 
and NGC 253 starbursts), the pionic secondary spectrum will also subside due to the kinetics of pion production. At low 
frequencies, such as those observed by LOFAR, both of these effects can increase the proton/electron ratio above that of the 
hadronic bremsstrahlung-losses case. At a few hundred MHz, where the secondaries are present, the ionization-loss formulae 
(section |2. 5. 2[) should work since it already takes into account the changing proton/electron ratio. The assumption of a CR 
proton injection spectrum of the form of a E~ p power law also probably starts to break down at proton energies of a GeV as 
well, with a momentum power law or some non-linear acceleration spectrum likely more appropriate. 

However, while these effects formally invalidate the bremsstrahlung-loss dominated formula outside of the frequency range 
where bremsstrahlung losses dominate (300 MHz to 6 GHz for M 82; 3 GHz to 15 GHz in Arp 220's nuclei), in practice, 
their effects are relatively benign. From the classical and revised formulae, the proton/electron ratio affects the magnetic 
field strength estimates only approximately as K 2 ^ 7 ' 2 : an order of magnitude change in Kq only results in a factor ~ 2 
change in _B m i n or B eq . At low frequencies, the ratio of ionization to bremsstrahlung loss times only goes as ti on /ibroms oc 
i/ 1 / 2 , and the proton/electron ra tio varies similarly. To take an extreme example, M 82's radio flux at 22.5 MHz is 39 Jy 
( Roger. Costain. fc Stewart 1986). Naively putting these values into the bremsstrahlung loss formula (and assuming all of 



this flux comes from the starburst region) gives _B cq = 170 fiG and B m i n — 160 fxG - values that are within ~ 25 per cent 
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of those obtained using the formula at 1 GHz. The more formally correct ionization-loss formulae instead give B cq = 180 fiG 
and B m i n = 140 fj,G. 



4 CONCLUSION 

We rederive minimum energy and equipartition estimates for the magnetic fields for the conditions that prevail in starburst 
galaxies. In these galaxies, the radio-emitting CR e population may consist largely of hadronic pp secondaries rather than 
primaries. Furthermore, strong radiative losses, particularly bremsstrahlung at GeV energies, set the lifetime of the CR e ± 
and affect the proton/electron ratio. Despite these effects, we found that the steady-state proton/electron ratio is probably 
close to ~ 100 for GHz-emitting e , assuming that ibrcms ~ 3t e . As a result, the classical and revised equipartition formulae 
give results that are quite similar to ours, despite their different assumptions. 

Although we scale the CR e ± lifetime to individual losses (particularly bremsstrahlung) to derive an analytical result, 
it should also be possible to numerically solve for _B cq and -B m i n using the full expression for i coo i including all losses, if 
equation [19] is used. However, since the density would not fully cancel out (due to synchrotron and IC losses in the 
lifetime), and because the IC lifetime depends on the radiation energy density, such an estimate would require both the gas 
density the CRs experience and the radiation energy density. 

While our focus has been on starburst galaxies, the formulae above should apply for any radio source where the CR e 



population is mostly pionic pp secondaries and is strongly cooled. iPfrom mcr & Enfilin ( 2004) already considered the case of 
radio emission from galaxy clusters, which may come from secondary e (|DennisorJll980f ) that cool from synchrotron and IC 
losses off the CMB over the Gyr life times of clusters. Active Galactic Nuclei may also accelerate CR protons when their jets 



are stalled by surrounding gas fe.g.. lAlvarez-Mufiiz fe Me szaros 2004). These protons may interact with the gas to produce 



pionic secondaries, which may radiate in the radio; the formulae derived here would also apply in such a case. 
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